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Fractal Risk Assessment of 1SS Propulsion Module in Meteoroid and 

Orbital Debris Environments 

Abstract 


A unique and innovative risk assessment of the International Space Station (ISS) 
Propulsion Module is conducted using fractal modeling of the Module s response to the 
meteoroid and orbital debris environments. Both the environment models and structural 
failure modes due to the resultant hypervelocity impact phenomenology, as well as 
Module geometry, are investigated for fractal applicability. The fractal risk assessment 
methodology could produce a greatly simplified alternative to current methodologies, such 
as BUMPER analyses, while maintaining or increasing the number of complex scenarios 
that can be assessed. As a minimum, this innovative fractal approach will provide an 
independent assessment of existing methodologies in a unique way. 

1.0 Introduction 


The meteoroid and orbital debris environments present considerable challenges to 
designers of large, long duration spacecraft, such as the ISS Propulsion Module This 
Module is a critical part of the ISS, performing altitude and attitude maintenance 
functions. Since it is a large piece of space hardware with a mission that lasts throughout 
ISS functionality, the Propulsion Module is expected to sustain and survive thousands of 
impacts by meteoroid and orbital debris particles. Therefore, a thorough assessment of its 
ability to withstand these environments is crucial to the success of the ISS mission overall 


Fractals are scale-invariant features that are becoming ubiquitous in the modeling of 
natural and man-made phenomena. Fractals are represented by simple mathematical 
models which display surprising degrees of complexity. In many cases, fractals have been 
found to more fully represent natural and engineering phenomena than much more 
complex physical and mathematical models Furthermore, fractals generally require much 
less data to “fit” the phenomena of interest. 

Typically, a meteoroid and orbital debris risk assessment is a complex undertaking 
Performing such an assessment generally requires expert information from a number of 
fields, including space environments, mission analysis, geometry modeling, structural 
design, and hypervelocity impact phenomenology. Assessments range from inexpensive 
“back of the envelope” calculations to moderately priced finite element approaches to 
expensive hydrodynamic analyses. In this field, there is always the question of how 
complex the assessment should be relative to the presumed and expected accuracy of the 
results There is also the issue of error propagation in the assessment. The more complex 
the assessment, the greater the potential for error at any step of the assessment Because 
the environments and structural response are highly nonlinear phenomena, any single error 
in the analysis can, in general, result in highly inaccurate results. 

Because the fractal approach results in greatly simplified model representations, both the 
costs to assess and the probability of assessment error are greatly reduced. Furthermore, 
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because simple fractal models result in complex phenomenology, the accuracy of the 
resulting assessment is not compromised, and could be enhanced. If fractals apply to this 
field, then a three-way improvement in meteoroid/orbital debris risk assessments could be 
realized At the very least, fractals will provide an alternative approach, not previously 
considered. 


The ISS Propulsion Module parameters used in this study are shown in Figure 1.0-1. 




i Orbital Parameters s 

Fixed Orbital Debris Particle 

- Average Altitude = 394.7 km 

Density = 2.8 gm/cm 3 * 

— Inclination = 51 .6 degrees & 

Effective Debris Area for 

- Operational from 10-1-2004 Through 

Module (L = 7 ft, D = 14 ft) = 

10-1-2014 

63.94 m 1 2 

, Whipple Shield Configuration 

- 0.063” Bumper (6061 -T6, 2.8 gm/cm 3 ) 

No Shadowing 

- 4” Bumper/Wall Standoff 


- 0.188” Wall (2219-T87, 62 ksi yield 

• '• O: . 

stress) 

// ,V\ \ ■- 

\ j \ ■ " 

Stuffed Whipple Shield Configuration 


- 0.075” Bumper (6061-T6, 2.8 gm/cm 3 ) 


- 4.5” Bumper/Wall Standoff 


- 0.188” Wall (221 9-T87, 62 ksi yield 


stress) 

- Standard Stuffing 

Tank Farm Minus Shielding 


Figure 1.0-1. ISS Propulsion Module Parameters for Meteoroid and Orbital Debris 

Risk Assessment. 


1.1 Study Objectives 

The objectives of this study are to investigate the applicability of fractal modeling to 
meteoroid/orbital debris risk assessment of spacecraft, and, if applicable, to demonstrate 
such using the ISS Propulsion Module The study will help answer the following key 
questions: 

1. Can fractal models accurately reflect complex relationships in the meteoroid and 

orbital debris environments 9 

A. Meteoroid (’racial is visually and statistical!} significant. 

B. Orbital debris fractal is visually and statistically significant. 
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2. Can fractal models accurately reflect complexities in geometry modeling of spacecraft 0 

A. Standard shapes are not fractals but are sell-similar. 

B. "Fractal" model has integer dimension. 

C'. Fractal "finite segment" analysis is meaningful at local level and identical in 
overall risk assessment at global level. 

3 Can fractal models accurately reflect complex phenomenology in hypervelocity 
impact? 

A. Whipple shield fractal is visually and statistically significant. 

B. Stuffed \\ hippie shield fractal is visually and statistically significant. 

4. How does a fractal risk assessment of the ISS Propulsion Module compare with an 
assessment using the traditional spreadsheet analysis? 

A. Relative error for orbital debris fractal risk assessment is *0.15% for 
Whipple Shield and -0.10% for Stuffed Whipple Shield. 

B. Relative error for meteoroid fractal risk assessment is for \N hippie 

Shield and -0.0014% for Stuffed Whipple Shield. 


1.2 Expected Significance 

The expected significance of this research and development effort includes an alternative 
approach to meteoroid/orbital debris risk assessments with the following potential 
advantages: 

1 . Reduced analysis complexity. 

2. Reduced analysis error, 

3. Increased understanding of risk and uncertainty. 

1.3 Statement of Uniqueness 

The fractal approach has not been previously applied to meteoroid/orbital debris risk 
assessments of spacecraft 

1.4 Current Schedule Status 


All items are completed. 


2.1 Task 1: Investigate Fractal Applicability to Meteoroid and Orbital 
Debris Environments 

Under this task, appropriate fractal models for the meteoroid and orbital debris 
environments are developed. The environments chosen are consistent with those to be 
experienced by the ISS Propulsion Module In this task, the environment models, orbital 
parameters, exposure times, and other applicable mission parameters are baselined relative 
to the ISS Propulsion Module. Recommendations are provided to MSFC personnel 
concerning all appropriate models and parameters that maintain consistency between 
fractal methodologies and BUMPER analyses 


2.1.1 Sample Meteoroid Environment Fractal Model 

As an example of the application of fractals to meteoroid/orbital debris environments, 
consider the portion of the meteoroid environment given by Cour-Palais [1969] as: 

Log l0 (N) = - 1 4.37 - 1 .2 1 3 Log w (m) 
me[l(T 6 .l] 

where, 

N = meteoroid flux (number of particles/m 2 /sec of mass m or greater), 
m = meteoroid particle mass (gm). 

This equation may be rewritten as: 

7V = 10~ 1437 m~ 1213 

This equation already satisfies the fractal power law property of scale invariance given by 


N = cm D 


where 

c = a constant, 

D = the fractal dimension. 

Thus, for this portion of the referenced meteoroid environment, the fractal dimension is 
1.213. This means that the dimension of the space of interest is between linear and 
quadratic, tending to be closer to linear While this fractal may or may not be appropriate 
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for use in the actual study, it is indicative of the presence of fractals in particulate space 
environments. 


The pertinent meteoroid environment applicable to the ISS Propulsion Module is given by 
Anderson and Smith [1994] The flux is given by: 



where 

F m (m) = meteoroid flux (number of particles/m 2 /year of mass m or greater), 
m = meteoroid mass (grams), 

Co = 3.156E07, 

Ci = 2.2E03, 
c 2 = 15, 
c, = I.3E-09, 

C4 = 1 0E1 1, 

c 5 = 1.0E27, 
c 6 = 1 3E-16, 

C 7 = 1 0E06. 


The density of the meteoroid particle varies with the mass of the particle. Under the 
assumption of a spherical particle, this density-mass variation corresponds to a density- 
particle density variation. The “break” points are given in Table 2 1 1 - 1 


Table 2.1. 1-1. Meteoroid Particle Density Break Points 


Particle Diameter (cm) 

Particle Mass (grams) 

Particle Density (gm/cm3) 

< 9.847E-03 

< 1.0E-06 

2.0 

> 9.847E-03 and < 0.2673 

>1.0E-06 and <0.01 

1.0 

> 0.2673 

>0.01 

0.5 


Fitting the above model to the simple fractal model, F m (m) - cm' D , we get: 


FJm) = (7.625£ - 08)w'' 17824 


The percentage variation in this model fit is given by R — 0.995, with significance level of 
1.79E-36, indicating an excellent statistical fit (Note that all R values reported are for the 
ln(Y) = ln(a) + b ln(X) model, not for the Y = a X b model.) Both estimated parameters 
also showed significance levels corresponding to an excellent fit. Thus, for the referenced 
meteoroid environment, the fractal dimension is approximately 1.178 This means that the 





dimension of the space of interest is between linear and quadratic, tending to be closer to 
linear 


The fit of this fractal to the actual environment definition is shown below in Figure 2.1.1- 

1 Notice that the fit is very good, except below particle diameters of 5.0E-03 cm, where 
the fractal tends to overestimate the flux significantly. Particle diameters in this range are 
not expected to do significant damage to the 1SS Propulsion Module, assuming pressure 
wall perforation is the critical failure mode of concern and that the module shielding is 
protective in this range of small particles. For the range of particle diameters between 1E- 

02 cm and 0.6 cm. the fractal model underestimates the meteoroid flux relative to the 
actual model, and above 0.6 cm, the fractal model tends to overestimate the flux relative 
to the actual model It is in this final range of particle diameters where the greatest 
damage potential lies for meteoroid impacts. Therefore, the fractal model would generally 
be considered to be conservative relative to the actual model. 
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Figure 2.1. 1-1. Comparison of Actual vs. Predicted (Fractal) Meteoroid 

Environment Model 


2.1.2 Orbital Debris Environment Fractal Model 

The pertinent orbital debris environment applicable to the ISS Propulsion Module is given 
by Anderson and Smith [1994] in NASA TM-4527. The equations for this model are 
given by: 
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F,(d,h,Ut,S ) = H(d)<Kh,S)'V(i)[F l (d)g l (t) + F,(rf)g 2 (r)] 

__ [J Qexp(-(log, 0 ^-0-78 )* 0.637 ) jl 2 

^(/2,5) = ^(^,5)/(^(/z,5) + 1) 

</> x (h,S) = 1 0 (/ ' 200-5 140-1 ' 5) 

F ] (d) = \.22x\0~ 5 d~ 25 

F 2 {d) = S.\x\0'°(d + 700y 6 

o+^r ,98 v<2oii 
g|( ' )_! (l + ?) 3 (l + g , r 20 ".'2 20 U 

g 2 (r) = l + p(;-1988) 


where, 

F r = flux, impacts per square meter of surface per year, 
d = orbital debris diameter (cm), 
t = date (year), 
h = altitude (km), 

S = 13-month smoothed solar radio flux for t-1 year (10 4 Jy), 
i = inclination (degrees), 

p = annual growth rate of mass in orbit (baseline = 0.05), 

q and q’ = estimated growth rates of fragment mass (0.02 and 0.04, respectively). 

For the ISS Propulsion Module, h = 394.7 km, i = 51.6 deg, and t = 10-1-2004 through 
10-1-2014 Figure 2.1 2-1 shows the resulting orbital debris flux for these parameters in 
the first year of exposure This figure also shows the corresponding flux from a computer 
model ORDEM96 for the same mission parameters. It should be noted that the 
ORDEM96 model appears to use 95%-ile solar flux values, while the flux shown for 
NASA TM-4527 is derived using 50%-ile solar flux values. Nevertheless, the flux values 
for the two environment definitions are comparable for particle diameters below about 
0.03 cm Above this threshold, however, the ORDEM96 values are significantly below 
those for NASA TM-4527 
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NASA TM-4527 
ORDEM96 


Diameter (cm) 


Figure 2. 1.2-1. Comparison of NASA TM-4527 vs. ORDEM96 Orbital Debris 
Environment Models. ESS Propulsion Module Parameteis (Orbital debris particles 
sweater than 10 cm are tracked and avoided.) 


Using the NASA TM-4527 orbital debris environment definition model, the total mission 
flux was calculated for the 1SS Propulsion Module mission and orbital parameters This 
was then fit with a simple fractal model of the form, 

Fn(d) 


cd' D , w e get: 


F d (d) = ( 2 . 8435 £ - 04)d 


- 2.1308 


The percentage variation in this model fit is given by R = 0.99, with significance level of 
] 5E-30, indicating an excellent statistical fit. (Note that all R values reported are for the 
ln(Y) = ’ln(a) + b ln(X) model, not for the Y = a X b model.) Both estimated parameters 
also showed significance levels corresponding to an excellent fit. Thus, for the referenced 
orbital debris environment, the fractal dimension is approximately 2.13. This means that 
the dimension of the space of interest is between quadratic and cubic, tending to be closer 
to quadratic. 

The fit of this fractal to the actual environment definition is shown below. Notice that the 
fit is very good, except above particle diameters of 4 cm, where the fractal tends to 
underestimate the flux significantly Particle diameters in this range may or may not 
significantly affect the risk calculations for the ISS Propulsion Module. If it is found that 
this range of particle diameters is a significant consideration, a multifractal (piecewise) 
approach will be adopted. In that case, the debris environment will be modeled by two or 
more fractals over various particle diameter ranges. For the range of particle diameters 
between IE-03 cm and IE-02 cm, the fractal model underestimates the debris flux relative 
to the actual model, and from IE-02 to 4 cm, the fractal model tends to overestimate the 
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flux relative to the actual model. It is in this final range of particle diameters where the 
greatest damage potential lies for orbital debris impacts. Therefore, the fractal model 
would generally be considered to be conservative relative to the actual model for particle 
diameters below about 4 cm 



(Orbital debris particles greater than 10 cm are tracked and avoided.) 


2.1.3 Graphical Display of Environment Fractal Models 


It was previously published that the orbital debris and meteoroid environment fractal 
dimensions pertinent to the ISS Propulsion Module are 1 1782 and 2.1308, respectively. 
How can we graphically display these dimensions in a 2-D way? One way is to consider 
the unit square: 
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We can replace each line (side) of this square by a series of lines smaller than the original. 
How do we determine the number of lines and the length of each 7 One way is to consider 
the equation: 


\_ 

r = N° 


where r = the similarity ratio (the factor by which the line exceeds the reduced lines), 

N = the number of “reduced” lines, and D = the dimension of the fractal 

But, how do we determine r and N, given only D 7 We desire the smallest values of r and 
N (to minimize complexity) that satisfy the equation above. Of course, there will always 
be some error associated with such a requirement for integer values Hopefully, we can 
find values that are relatively small, yet produce small errors. In fact, this was achieved 
for these environments. 


For the meteoroid environment (D — 1.1782), we find associated values of r 4 and N 
5, with an error of only 2% in the determination of r This means that we can replace each 
line of the square with 5 lines, each 1/r = V4 the original length. There are a number of 
ways to achieve this, but we would like to preserve as much symmetry as possible One 
method of replacement (formally called a generator) is shown: 
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Note the symmetry about a vertical line for this replacement scheme. Following the initial 
replacement, each successive (reduced) line is replaced by 5 lines, each V* the original 
length. This process is continued infinitely many times, or until further resolution is 
impossible. Figure 2. 1 .3-1 shows the resulting 2-D fractal graphic. This graphic shows 
four dense regions linked by less dense regions. Notice that the hole in each link 
represents an identical but much smaller representation of the interior of this graphic 
Many other self-similarities are visible in this figure. What technical conclusions, it any, 
can we draw from such graphical depictions of the meteoroid fractal environment for the 
1SS Propulsion Module 9 


For the orbital debris environment (D — 2 1308), we find associated values of r 2 and N 
= 4, with an error of only 4% in the determination of r. This means that we can replace 
each line of the square with 4 lines, each 1/r = 1/2 the original length. Again, there are a 
number of ways to achieve this, but we would like to preserve symmetry if possible A 
three dimensional representation is preferred for a fractal with dimension greater than 2. 
However, we can get some graphical feel for this orbital debris fractal in two dimensions. 
One method of replacement for this fractal is shown: 
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Again, note the symmetry about a vertical line for this replacement scheme Following the 
initial replacement, each successive (reduced) line is replaced by 4 lines, each 1/2 the 
original length This process is continued infinitely many times, or until further resolution 
is impossible. Figure 2.1 3-2 shows the resulting 2-D fractal graphic for !4 of the total 
graphic. This graphic shows much higher interior density than that for the meteoroid 
fractal This is indicative of the higher dimensionality of the orbital debris environment 
fractal, relative to that for the meteoroid environment. Notice the cascading identical 
interiors in the “claw” patterns of this graphic. Many other self-similarities are visible in 
this figure. What technical conclusions, if any, can we draw from such graphical 
depictions of the orbital debris fractal environment for the 1SS Propulsion Module 9 
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Figure 2. 1.3-1. "Smilm.i, ! rr/“\ lioltiui- Homh" Fractal Representation of 
Meteoroid Environment Associated With 1SS Propulsion Module (Dimension — 
1.1782). Image created by replacing each line of a square by 5 lines, each Vi the length of 
the original, according to the formula r = (N) 1 D where r = scaling factor (4), N = number 
of lines replacing original line (5), and D = dimension of fractal (1.1 782). Error - 2% 
(Image created using BRAZIL DESIGN software.) 
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Figure 2.1.3-2. "Debris Clou's" One-Fourth of 2-D Fractal Representation of 
Orbital Debris Environment Associated With ISS Propulsion Module at Altitude of 
394.7 km (Dimension = 2.1308). Image created by replacing each line of a square by 4 
lines, each 1/2 the length of the original, according to the formula r = (N) where r = 
scaling factor (2), N = number of lines replacing original line (4), and D = dimension of 
fractal (2.1 308). Error = 4%. Note the higher “interior” density of the orbital debris 
fractal, as compared to the meteoroid environment fractal. (Image created using BRAZIL 

DESIGN software.) 



2.2 Task 2: Investigate Fractal Applicability to Geometry Modeling 

In this task, appropriate fractal models for geometry modeling of spacecraft are 
developed. The models chosen are consistent with the ISS Propulsion Module 
Information on the ISS Propulsion Module geometry, bumper/wall configuration, and 
BUMPER geometry model is collected and reviewed. Recommendations are provided to 
MSFC personnel concerning all appropriate models and parameters that maintain 
consistency between fractal methodologies and BUMPER analyses. 


2.2.1 Sample Geometry Fractal Model 

As a simple example of the recursive, self-replicating application of fractals to geometry 
modeling, consider the simple rectangle below. The rectangle is representative of a side 
view of a cylindrical shaped module, where L is the length of the cylinder and D is the 
diameter (generally, external) of the cylinder. 




4 

-*| 


D 


L 


Figure 2.2.1-1. Oth Stage of Fractal Geometry Model 

This rectangle might represent, for instance, a highly simplified cross-sectional view of the 
ISS Propulsion Module. If the recursion rule for this geometry is to halve the rectangle(s) 
at each step, the second step geometry would look like: 
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◄ 

► 


L 

Figure 2.2. 1-2. First Stage of Fractal Geometry Model 





At the fourth stage, the geometry would appear as follows: 



Figure 2.2.1-3. Fourth Stage of Fractal Geometry Model 

This stage of the fractal may better represent an idealized version of the cross-section of 
the ISS Propulsion Module. Other, more complicated examples of geometry modeling 
can be expected 

Note that the number of subdivided areas (rectangles), N(A n ), is given by: 


16 




N(A h ) = 2 " , m = 0 , 1 , 2 , 3 ,... 


where n is the stage of subdivision, and A n is the area of each individual subarea at this 
stage The case, n = 0, corresponds to the single original rectangle of area LD Clearly, 
there is one such rectangle of area LD, as indicated by the formula. The case, n 1, 
corresponds to 2 rectangles and the case n = 4 corresponds to 16 rectangles, as given by 
the formula 

Furthermore, the area of each subdivided rectangle at the nth stage. A„, is given b\ 



The case n = 0 corresponds to the single original rectangle of area LD, n - 1 corresponds 
to two rectangles each of area LD/2, and the case n = 4 corresponds to 1 6 rectangles each 
of area LD/16, all correctly predicted by the formula. 

Several interesting properties emerge from these simple formula. The first is that 



Thus, regardless of the number of subdivisions, the product of the number of subareas and 
the area of each individual subarea is equal to the total original area, LD This equation 
may be rewritten as: 

N(A„) = LD(A„T',n= 0 , 1 , 23 ,-. 


This fits the fractal form with c = LD, and dimension - 1. Clearly, this is not a fractional 
exponent, so the term “fractal” is generalized (to include integer values). Nevertheless, 
the self-similarity of the rectangular subareas is consistent with fractal use 

Also, by the property of a geometric series, the infinite sum of these subareas is given by: 



Thus, it follows that the partial sum is given by: 
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±A„=± a „~ a 0 = 2 LD-^ = 2LD-LD= LD 

=0 


In this sense, the “fractal” use of self-similar geometric objects may be viewed as either an 
infinite sum or the product of the number of subareas and the area of each subarea at a 
given stage of subdivision Spacecraft shadowing can also be taken into account b\ 
limiting the number of areas, N(A„), at a given stage of subdivision. 


A schematic of the ISS Propulsion Module is shown in Figure 2.2. 1-4 Projected areas of 
this configuration can be dealt with in three ways: the standard BUMPER analysis 

method, or the fractal method, or some combination of both Each method will be 
investigated during the geometry, environment, and phenomenology tasks tor the 
meteoroid and/or orbital debris environments 




Figure 2,2. 1-4. Sketch of ISS Propulsion Module - Tank Farm Minus Shielding 

(Courtesy of NASA-MSFC) 
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2.3 Task 3: Investigate Fractal Applicability to Hypervelocity Impact 
Phenomenology Modeling 


This section involves the development of appropriate fractal models for hypervelocity 
impact phenomenology modeling of spacecraft. The models chosen are consistent with 
those applicable to the ISS Propulsion Module. Information on meteoroid/orbital debris 
critical failure modes, existing penetration equations, and other appropriate quantitative 
phenomenological models (eg., crater area, depth) are reviewed and chosen for 
applicability and comparison to those provided in BUMPER. Recommendations are 
provided to MSFC personnel concerning all appropriate models and parameters in order 
to maintain consistency between fractal methodologies, spreadsheet approaches, and 
BUMPER analyses. 


2.3.1 Whipple Shield Phenomenology Fractal Model (Orbital Debris) 

A fractal for a Whipple Shield phenomenology model relevant to ISS Propulsion Module 
is developed in this section 

The Whipple Shield Ballistic Limit model (Christiansen, 1993) is given over three impact 
velocity ranges as: 


d c =[(C(^)° 5 +O/(O.6(cos0) 5/ V,°, 5 F 2 ' 3 )f /l 
V cosO < 3 km / see 


d =! [ ( , u . ( ^)» 5 +^)/(1.248p; 5 cos6>)f' ,l '(l.75-(l'cos0)/4)! 

+ {[i.07U; , V; l 'V; l,, s l ' 5 (^) 1 ' 3 ]((r , cos^)/4-0 75)} 

3 km / sec < V cos 6 < 7 km / sec 
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id c =3.918^V; l/3 A 1/9 (Fcos0r 2/3 5'' 3 (^) ,/3 


V cos 6 > 1km / sec 


where: 

d c = critical projectile diameter (cm) causing failure, 
p = density (g/cirf ), 

S = overall spacing between outer bumper and rear wall (cm), 
o = rear wall yield stress (ksi), 
t = thickness (cm). 

0 = impact angle (deg) measured from surface normal, 

V = projectile velocity (km/sec), 

with subscripts: 
b = bumper, 
p = projectile, 
w = rear wall. 

The configuration analyzed is a Whipple Shield with 0.063” bumper, 4 spacing, and 
0.188” rear wall. A ballistic limit curve (critical particle diameter vs. impact velocity) is 
developed for the orbital debris environment with fixed particle density of 2.8 gm/cm (as 
specified in the ISS design requirements) using the equations above. A fractal is then fit to 
this ballistic limit equation The best fit fractal is given by: 



2.917 V ~ 0 6028 


This fractal shows an excellent fit to the actual ballistic limit curve. The correlation (R) of 
0.991 together with significance level of 4E-54 indicates a remarkably good fit by this 
fractal (Note that all R values reported are for the ln(Y) = ln(a) + b ln(X) model, not for 
the Y = a X b model ) Furthermore, the visual fit shown below in Figure 2.3. 1-1 is 
excellent. 
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Figure 2.3.1-1. Comparison of actual ballistic limit curve for Whipple shield over 
orbital debris velocity range with derived fractal model. Both statistical and visual 
significance indicate an excellent fractal fit. 1SS Propulsion Module Parameters 


Figure 2.3. 1-2 shows the graphical representation for the Whipple Shield phenomenology 
fractal: 


D cnl = 2.9 17F“° 6028 


This graphic was obtained by replacing each line in the unit square by 4 lines, each 1/10 
the length of the original line. The relative sparseness of this representation, as compared 
to the meteoroid and orbital debris flux fractals is indicative of the low fractal dimension, 
0.6028 


Actual 

Fractal 
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Figure 2.3.1-2. / / .. . r\:urn’ 2-D Fractal Representation of Whipple Shield 

Phenomenology for Orbital Debris Associated With ISS Propulsion Module at 
Altitude of 394.7 km (Dimension = 0.6028). Image created by replacing each line of a 
unit square by 4 lines, each 1/10 the length of the original, according to the formula r = 
(N) in where r = scaling factor (10). N = number of lines replacing original line (4), and D 
= dimension of fractal (0.6028). Error = 0.3%. Note the lower “interior” density of the 
phenomenology fractal, as compared to the meteoroid and orbital debris flux fractals. 

(Image created using BRAZIL DESIGN software.) 
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2.3.2 Stuffed Whipple Shield Phenomenology Fractal Model (Orbital Debris) 

The equations for the Stuffed Whipple Shield configuration are given in Christiansen et al 
[1995], The fractal for the Stuffed Whipple Shield configuration (orbital debris 
environment) of the ISS Propulsion Module is given by: 



86.73K” 1 789 


Figure 2. 3. 2-1 shows the visual fit of this fractal (predicted) to the actual data values. The 
correlation (R) of this fractal model is 0 9766, with a significance level of 9 5E-42. 
indicating a statistically excellent fit (Note that all R values reported are for the ln(Y) = 
ln(a) + b ln(X) model, not for the Y = a X b model.) The visual fit also appears to be quite 

good. 

Figure 2.3. 2-2 shows the graphical representation (1/4 only) for the Stuffed Whipple 
Shield phenomenology fractal This graphic was obtained by replacing each line in the unit 
square by 7 lines, each 1/3 the length of the original line. 
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Figure 2.3.2-1. Phenomenology Fractal for Stuffed Whipple Shield Configuration of 
ISS Propulsion Module. ISS Propulsion Module Parameters 





Figure 2.3.2-2. j'-aiilin v, C n> su-.s 2-D Fractal Representation of Stuffed Whipple 
Shield Phenomenology for Orbital Debris Associated With ISS Propulsion Module 
at Altitude of 394.7 km (Dimension = 1.789). Image created by replacing each line of a 
unit square by 7 lines, each 1/3 the length of the original, according to the formula r = 
(N) id where r = scaling factor (3), N = number of lines replacing original line (7), and D = 
dimension of fractal ( 1 789) Error = 1%. (Image created using BRAZIL DESIGN 

software.) 
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2.3.3 Sample Crater Phenomenology Fractal Model 

Klop et al [1990] present a graph of the fraction of witness plate craters following 
hypervelocity impact versus the crater diameter Approximate values from this graph are 
shown below in Table 2. 3. 3-1. 


Table 2.3.3-1 Approximate Values from Klop et al. [1990] for Experiment 11. 


Crater Diameter (mm) 

Fraction of Craters With Crater Diameter ! 

0.2 

0.15 ! 

0.3 

0.205 

0.4 

0.14 

0.5 ^ 

0.1 

0.6 

0.085 

0.7 

0.08 

0.8 

0.06 

0.9 1 

0.035 

1.0 

0.03 

1.1 

0.025 

1.2 

0.022 

1.3 

0.015 

1.4 

0.01 

1.5 

0.004 

1.6 

0.02 

1.7 

0.004 

1.8 

0.0001 

1.9 

0.004 

2.0 

0.0001 

2.1 

" 0.004 

2.2 

0.0018 


The power law for fractals in this case is given by: 



where 

f c = the fraction of craters with diameter at least as great as d c , 
a = a constant, 

dc = the crater diameter (mm), 

D = the fractal dimension 

Regression of Table 2 3 3-1 values results in the following fractal equation: 
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f c = 0.07 Ud; 2 505 
d c e [0.35,2.25] 


The correlation for this regression is R = 0.922, while the significance level is a ’ 03E- 
09 indicating a highly sienificant model from a statistical standpoint (Note that all R 
values reported are for the ln(Y) = ln(a) + b ln(X) model, not for the Y - a X model ) 
Using this model, we see that the fractal dimension for crater diameter features on the 
witness plate is approximately 2.5. This means that the dimension of the space ot crater 
diameters is roughly halfway between quadratic (2-D) and cubic (3-D). Note, too, that if 
the equation above is changed to model the total number of craters instead of the fraction 
of craters, the coefficient will change, but the exponent will not Thus, the fractal 
dimension will remain the same regardless of the method of data normalization Figure 
2. 3. 3-1 shows this relationship The area of the bubbles, while not to scale, represent the 
relative projected areas of the craters for a given crater diameter While this fractal may 
or may not be appropriate for use, it is indicative of the presence of fractals in 
hypervelocity impact phenomenology models other than those for ballistic limit curves 

Figure 2. 3. 3-2 shows the corresponding graphical representation for Vi of this fractal 
Note the relatively higher fractal density of this graphic. This is due to the high dimension 
(2.505) associated with this fractal 



Figure 2.3.3-1. Fraction of Craters Greater Than Given Crater Diameters 
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Figure 2. 3. 3-2. "I thm'ii Smiling r ro^ ! n\\ u/" Hiccns 2-1) Fractal 

Representation of Crater Phenomenology (Dimension = 2.505). Image created by 
replacing each line of a unit square by 6 lines, each 1/2 the length of the original, 
according to the formula r = (N) n> where r = scaling factor (2), N = number of lines 
replacing original line (6), and D = dimension of fractal (2.505). Error — 2 o. Greater 
fractal density represents higher fractal dimensionality. (Image created using BRAZIL 

DESIGN software.) 


2.3.4 Whipple Shield Phenomenology Fractal Model (Meteoroid Environment) 

A fractal for a Whipple Shield phenomenology model relevant to 1SS Propulsion Module 
in the meteoroid environment is developed in this section. 

The Whipple Shield Ballistic Limit model (Christiansen, 1993) is given over three impact 
velocity ranges as: 




d c = [(* 

V ms ft 


(£>■ 5 + t h ) /(0.6(cos df : \p° p 5 V 213 )] 
< 3km / sec 


18/19 


d e ={[(t w (^) 05 +t h )/(\.2ASp° p 5 cos#)] 18/l9 (1.75 - (V cos^)/4)} 


+ {[\.Ql\t™ cos 6>)/4 -0.75)} 

3 km / sec < V cos 9 < 1km /sec 


lJ c =3.918r„ 2/ V; 1/3 A 1,9 (Fcos0r 2/3 5 1/3 (^) 1/3 


V cos 6 > 1km / sec 


where: 

d c = critical projectile diameter (cm) causing failure, 
p = density (g/cm'), 

S = overall spacing between outer bumper and rear wall (cm), 
a = rear wall yield stress (ksi), 
t = thickness (cm). 

0 - impact angle (deg) measured from surface normal, 

V = projectile velocity (km/sec), 

with subscripts: 
b = bumper, 
p = projectile, 
w = rear wall. 

The configuration analyzed is a Whipple Shield with 0.063 bumper, 4 spacing, and 
0.188” rear wall. A ballistic limit curve (critical particle diameter vs impact velocity) is 
developed for the meteoroid environment using the equations above A fractal is then fit 
to this ballistic limit equation The best fit fractal is given by: 


28 





d„ = 7.016 V 


- 0.8143 


( 0 + 1 ) 


0.1595 


This fractal shows an excellent fit to the actual ballistic limit curve. The correlation (R) of 
0.886 together with significance level of 3E-144 indicates a remarkably good fit by this 
fractal (Note that all R values reported are for the ln(Y) = ln(a) + b ln(X) model, not for 
the Y = a X b model.) The visual fits for various values of impact angle are shown below 

in Figures 2.3 .4-1 through 2. 3. 4-7. 



Figure 2.3.4- 1. Comparison of actual ballistic limit curve for Whipple shield over 
meteoroid velocity range (0 degree impact angle) with derived fractal model. 

* ISS Propulsion Module Parameters 
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Figure 2.3.4-2. Comparison of actual ballistic limit curve for Whipple shield over 
meteoroid velocity range (15 degree impact angle) with derived fractal model. 

!SS Propulsion Module Parameters 



Figure 2.3.4-3. Comparison of actual ballistic limit curve for Whipple shield over 
meteoroid velocity range (30 degree impact angle) with derived fractal model. 

ISS Propulsion Module Parameters 
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Figure 2.3.4-4. Comparison of actual ballistic limit curve for Whipple shield over 
meteoroid velocity range (45 degree impact angle) with derived fractal model. 
ISS Propulsion Module Parameters 




Figure 2.3.4-5. Comparison of actual ballistic limit curve for Whipple shield over 
meteoroid velocity range (60 degree impact angle) with derived fractal model. 

ISS Propulsion Module Parameters 
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Figure 2.3.4-6. Comparison of actual ballistic limit curve for Whipple shield 
over meteoroid velocity range (75 degree impact angle) with derived fractal 
model. ISS Propulsion Module Parameters 



Figure 2.3.4-7. Comparison of actual ballistic limit curve for Whipple shield over 
meteoroid velocity range (90 degree impact angle) with derived fractal model. 

ISS Propulsion Module Parameters 


Figure 2. 3. 4-8 shows the graphical representation for the Whipple Shield phenomenology 
fractal: 
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A™ =7.0161' 


-0.8143 


(0 + 1 ) 


0.1595 


This graphic was obtained by replacing each line in the unit square by 6 lines, each 1/9 
the length of the original line The fractal variable shown is velocity, with fractal 
dimension of 0 8143 


Figure 2.3.4-8. in ic of t irclcs 2-D Fractal Representation of Whipple Shield 

Phenomenology for Meteoroids Associated With ISS Propulsion Module at Altitude 
of 394.7 km (Dimension = 0.8143). Image created by replacing each line of a unit square 
by 6 lines, each 1/9 the length of the original, according to the formula r = (N) where r 
= scaling factor (9), N = number of lines replacing original line (6), and D = dimension of 
fractal (0.8143). Error = 0.33%. (Image created using BRAZIL DESIGN software.) 




2.3.5 Stuffed Whipple Shield Phenomenology Fractal Model (Meteoroid 
Environment) 

The equations for the Stuffed Whipple Shield configuration are given in Christiansen et al 
[1995], The fractal for the Stuffed Whipple Shield configuration (meteoroid environment) 
of the ISS Propulsion Module is given by: 


D cril = 6.8 W -° 5319 (6> + l)° 2679_ 


Figures 2.3 5-1 through 23.5-7 show the visual fits of this fractal to the actual data values 
for various particle impact angles. The correlation (R) of this fractal model is 0.515, with 
a significance level of 1.4E-29, indicating a statistically excellent fit. (Note that all R 
values reported are for the ln(Y) = ln(a) + b ln(X) model, not for the Y = a X model ) 
The main reason for the relatively low R value is the exceedingly high critical particle 
diameters for low velocities and impact angle of 89 degrees. 

Figure 23.5-8 shows the graphical representation for the Stuffed Whipple Shield 
phenomenology fractal. This graphic was obtained by replacing each line in the unit 
square by 3 lines, each 1/8 the length of the original line. 
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Figure 2.3.5-1. Phenomenology Fractal for Stuffed Whipple Shield Configuration of 
ISS Propulsion Module in Meteoroid Environment (0 degree impact angle). 
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Figure 2.3.S-2. Phenomenology Fractal for Stuffed Whipple Shield Configuration of 
ISS Propulsion Module in Meteoroid Environment (15 degree impact angle). 
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Figure 2.3.5-3. Phenomenology Fractal for Stuffed Whipple Shield Configuration of 
ISS Propulsion Module in Meteoroid Environment (30 degree impact angle). 

ISS Propulsion Module Parameters 


45 Degree Impact Angle 


5.000 



0.000 


-♦—Actual 
■t Fractal 


0 20 40 60 80 

Particle Velocity (km/sec) 


100 


Figure 2.3.5-4. Phenomenology Fractal for Stuffed Whipple Shield Configuration of 
ISS Propulsion Module in Meteoroid Environment (45 degree impact angle). 
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Figure 2.3.5-5. Phenomenology Fractal for Stuffed Whipple Shield Configuration of 
ISS Propulsion Module in Meteoroid Environment (60 degree impact angle). 

ISS Propulsion Module Parameters 
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Figure 2.3.S-6. Phenomenology Fractal for Stuffed Whipple Shield Configuration of 
ISS Propulsion Module in Meteoroid Environment (75 degree impact angle). 
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Figure 2.3.5-7. Phenomenology Fractal for Stuffed Whipple Shield Configuration of 
ISS Propulsion Module in Meteoroid Environment (90 degree impact angle). 

ISS Propulsion Module Parameters 
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Figure 2.3.S-8. "! nawnii -Jr;!'" 2-D Fractal Representation of Stuffed W'hipple 
Shield Phenomenology’ for Meteoroids Associated With ISS Propulsion Module at 
Altitude of 394.7 km (Dimension = 0.5319). Image created by replacing each line of a 
unit square by 3 lines, each 1/8 the length of the original, according to the formula r = 
(N) 1 D where r = scaling factor (8), N = number of lines replacing original line (3). and D = 
dimension of fractal (0.53 1 9) Error = 1 .4%. (Image created using BRAZIL DESIGN 

software.) 



2.4 Task 4: Comparison of Fractal Risk Assessment With BUMPER Results 

The purpose of this task is to compare and contrast the fractal risk assessment of the ISS 
Propulsion Module with that provided by the BUMPER code and traditional spreadsheet 
risk assessment. The spreadsheet inputs and outputs are compared for consistency with 
the fractal approach The fractal risk assessment is executed with consistency and 
completeness in mind. The results of the two methods (spreadsheet and fractal) are 
compared and contrasted statistically and qualitatively. Recommendations are made 
regarding improvements in the fractal approach A “finite segment’' analysis is performed 
using both spreadsheet and fractal methods, and producing a “risk plot.” showing the risk 
to orbital debris particles on the front face and sides of the Propulsion Module This 
preliminary “risk plot” may be useful in comparison with equivalent BUMPER results. 


2.4.1 Orbital Debris Risk Assessment - Whipple Shield Fractal Model 


In this section, a risk assessment is performed for the ISS Propulsion Module with 
Whipple Shield configuration subject to the orbital debris environment with fixed particle 
density of 2.8 am/cm’. Figure 2.4. 1-1 shows the fit for the following local orbital debris 
flux fractal: 


F D (d) = 1 . 0509 E - 04 


This fractal was developed by considering the reduced range of values for critical orbital 
debris particle diameter pertinent to the Whipple Shield configuration of the ISS 
Propulsion Module By considering only these particle diameters, a more accurate “local 
orbital debris flux fractal can be obtained for analysis purposes. The correlation (R) for 
this fractal is 0.967, with significance level of 5.6E-36, indicating an excellent statistical fit. 
(Note that all R values reported are for the ln(Y) = ln(a) + b ln(X) model, not for the Y = 
a X b model.) 

The visual fit, while somewhat satisfactory could be improved by considering two ranges 
of particle diameters: above and below 1 cm, respectively. Nevertheless, we will see that 
the use of this fractal results in a high precision estimate of PNP, when compared to the 
equivalent traditional assessment. 

Figure 2.4. 1-2 shows the graphical representation of the local debris flux fractal. This 
fractal is obtained by replacing each line of a unit square by 7 lines, each !4 the length of 
the original line. The figure shows only % of the traditional “coastline” fractal for the 
debris flux If, instead of replacing each line of the square, we replace the entire square by 
7 squares, each l A the size of the original, we obtain the “gasket” fractal shown in Figure 
2.4. 1-3. 
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Figure 2 4 1-4 shows the fit for the following “non-fractal” approximating model to the 
orbital debris velocity probability distribution: 


f'(V) = 8.177£-04(K 14172 ) 


While the correlation (R) to the actual model of 0.985 is statistically significant (1 3E-45), 
the visual fit for this model leaves something to be desired. Nevertheless, we will see that 
employing this model in the integrated analysis is satisfactory for this particular situation 


Integration: 

We first combine the Whipple Shield Phenomenology fractal 

D cnl = 2.91 1V ~° 6028 

with the local orbital debris flux fractal 

FJd) = 1 ,0509£ - 04(iT 14164 ) 


to obtain the following (setting d = D cnt ): 

F„(V) = 1.0509 F- 04[2.917(K' , ’“ 2 * )]■' ■"“ 
= 2.3069£-05(l /M5,1 ‘) 


Next, we multiply this by the approximating model for orbital debris velocity probability 
distribution: 


f\V)F n (V) - 8. 1 11 E - 04CI 7 ' 4,72 )2.3069£ - 05(F° 8538 ) 
= 1.8863£’- 08(K 2 271 ) 


Now, we can calculate the ISS Propulsion Module orbital debris PNP for an effective area 
of 63.94 m 2 as: 
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15 


! PNPom, =e 


-63.94(4.244) Jl.8863£-08(F 2:!71 )dV 


3.271 


= e 


-5.1 1 87£-06(— ) 

3 271 = £-i IE-02 = 0.9891 


Note that the factor of 4.244 in front of the integral is necessary to guarantee that 


15 

4.244 jf'(V)dV = 1 

0 


The upper limit of both integrals corresponds to the highest (non-zero) velocity value for 
orbital debris particles for the situation analyzed. The equivalent orbital debris PNP using 
traditional risk assessment is 0.9876, resulting in a relative difference of 0. 1 5% 



Figure 2.4. 1-1. Local Orbital Debris Flux Fractal for Relevant Whipple Shield 
Particle Diameters. Fractal Model is “Predicted ” ISS Propulsion Module Parameters 
(Orbital debris particles greater than 10 cm are tracked and avoided.) 
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Figure 2. 4.1-2. • Liuhnwi"" 2-D Fractal Representation of Whipple 

Shield Pertinent Orbital Debris Flux for 1SS Propulsion Module at Altitude of 394.7 
km (Dimension = 1.4164). Image created by replacing each line of a unit square by 7 
lines, each 1/4 the length of the original, according to the formula r = (N) 1 D where r = 
scaling factor (4), N = number of lines replacing original line (7), and D = dimension of 
fractal ( 1 .4164) Error = 1 .3%. (Image created using BRAZIL DESIGN software.) 
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Figure 2.4.1-3. .“7SV% //;»" 2-D “Gasket” Fractal Representation of Whipple 

Shield Pertinent Orbital Debris Flux for ISS Propulsion Module at Altitude of 394.7 
km (Dimension = 1.4164). Image created by replacing each unit square by 7 squares, 
each 1/4 the size of the original, according to the formula r = (N) 1 D where r = scaling 
factor (4), N = number of squares replacing original square (7), and D = dimension of 
fractal ( 1 4164) Error = 1 .3%. (Image created using BRAZIL DESIGN software.) 
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2.4.2 Orbital Debris Risk Assessment - Stuffed Whipple Shield Fractal Model 

In this section, a risk assessment is performed for the ISS Propulsion Module with Stuffed 
Whipple Shield configuration subject to the orbital debris environment with fixed particle 
density of 2.8 gm/cm\ Figure 2.4 .2-1 shows the fit for the following “local” orbital debris 
flux fractal 


F„(d) = 8.4784E-050T' 0626 ) 


This fractal was developed by considering the reduced range of values for orbital debris 
critical particle diameter pertinent to the Stuffed Whipple Shield configuration of the ISS 
Propulsion Module By considering only these particle diameters, a more accurate “local” 
orbital debris flux fractal can be obtained for analysis purposes. The correlation (R) for 
this fractal is 0.931, with significance level of 6.8E-28, indicating an excellent statistical fit 
(Note that all R values reported are for the In(Y) = ln(a) + b ln(X) model, not for the Y = 
a X b model.) 

The visual fit, while somewhat satisfactory could be improved by considering two ranges 
of particle diameters: above and below 300 cm, respectively. Nevertheless, we will see 
that the use of this fractal results in a high precision estimate of PNP, when compared to 
the equivalent traditional assessment 
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Figure 2 4.2-2 shows the graphical representation of the local debris flux fractal This 
fractal is obtained by replacing each line of a unit square by 8 lines, each 1/7 the length of 
the original line. The figure shows only 14 of the traditional “coastline” fractal for the 
debris flux The relative smoothness of this fractal is indicative of a dimension in the 
vicinity of 1.0 (actual fractal dimension = 1 0626) 

Integration: 


We first combine the Stuffed Whipple Shield Phenomenology fractal 



86.73K” 1 789 


with the local orbital debris flux fractal 

F„(d) = 8.4784E - 05frr'° 626 ) 


to obtain the following (setting d = D cn t): 

F n (V) = 8.4784 £ - 05[86.73(K- 1789 )] I W,2< ’ 
= 7.3929 £-07(K 1901 ) 


Next, we multiply this by the approximating model for orbital debris velocity probability 
distribution: 


f\V)F D {V) = 8. 1 HE - 04(K 14,72 )7.3929/: - 07(F' 901 ) 
= 6.0452A;- 10(K 3 3,8 ) 


Now, we can calculate the 1SS Propulsion Module orbital debris PNP for an effective area 
of 63.94 m 2 as: 


15 


PNP l)ehn , = e 


-63.94( 4.244) J 6.0452/;-! 0( V y ~ ] *)dV 


15 


4.318 


) 


-1.6404/:-07(- , M ^ ^ _ 

= £> 4.318 =e -4.55E-O3 =09955 


Note that the factor of 4.244 in front of the integral is necessary to guarantee that: 
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15 

4 . 244 \f'(V)dV = 1 

0 


The upper limit of both integrals corresponds to the highest (non-zero) velocity value for 
orbital debris particles for the situation analyzed. The equivalent orbital debris PNP using 
traditional risk assessment is 0.9965, resulting in a relative difference of 0. 10 /o 



Figure 2.4.2-1. Local Orbital Debris Flux Fractal for Relevant Stuffed Whipple 
Shield Particle Diameters. Fractal Model is “Predicted 1S> Propulsion Module 
Parameters (Orbital debris particles greater than 10 cm are tracked and avoided ) 
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Figure 2.4.2-2. "u<‘iiinr HUh 2-D Fractal Representation of Stuffed Whipple 
Shield Pertinent Orbital Debris Flux for ISS Propulsion Module at Altitude of 394.7 
km (Dimension = 1.0626). Image created by replacing each line of a unit square by 8 
lines, each 1/7 the length of the original, according to the formula r = (N) 1 11 where r = 
scaling factor (7), N = number of lines replacing original line (8), and D = dimension of 
fractal (1.0626). Error = 1 . l° o (Image created using BRAZIL DESIGN software ) 


2.4.3 Meteoroid Risk Assessment - Whipple Shield Fractal Model 


In this section, a risk assessment is performed for the ISS Propulsion Module with 
Whipple Shield configuration subject to the meteoroid environment. Figure 2 4.3-1 shows 
the fit for the following “local” meteoroid flux fractal: 


F kt (d) = 6.357£-06(t/“ 3 333 )j 




This fractal was developed by considering the reduced range of values for critical 
meteoroid particle diameter pertinent to the Whipple Shield configuration of the 1SS 
Propulsion Module. By considering only these particle diameters, a more accurate “local’ 
meteoroid flux fractal can be obtained for analysis purposes. The correlation (R) for this 
fractal is 0.985, with significance level of <lE-235, indicating an excellent statistical fit 
(Note that all R values reported are for the ln(Y) = ln(a) + b ln(X) model, not for the T - 
a X b model.) Note that the impact angle is not a significant predictor variable for this 

model. 

Figure 2. 4. 3-2 shows the graphical representation of the local meteoroid flux fractal This 
fractal is obtained by replacing each line of a unit square by 10 lines, each 1/2 the length of 
the original line The figure shows only '/4 of the traditional “coastline” fractal for the 
meteoroid flux. It is very difficult to obtain a meaningful 2-D representation of a great er- 
than-3-D fractal 

Figure 2.4. 3-3 shows the fit for the following fractal model to the meteoroid velocity 
probability distribution: 


f'(V) = 2.653£-06(l / ' 556 ) 


The correlation (R) to the actual model of 0.986 is statistically significant (<lE-235), and 
the visual fit for this model is quite good Furthermore, we will see that employing this 
model in the integrated analysis is satisfactory for this particular situation. No meaningful 
fractal 2-D graphical representation can be made for a fractal with dimension greater than 

5. 


Integration: 

We first combine the Whipple Shield Phenomenology fractal 

D,„, = 7.016K 08l43 (6» + 1) 015,>5 

with the local meteoroid flux fractal 

~F m (d) = 6.35 7£ - 06(flf~ 3333 ) 

to obtain the following (setting d = D cn i): 

F u (Vj) = 6.357 £ - 06[7.016(F~° 8143 )(0 + if 595 ]' 3 ” 3 ' 
= 9.62 E - 09(V 2 7,4 ){0 + l)-° 53,6 
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Next, we multiply this by the fractal model for meteoroid velocity probability distribution 

7(F)F u {V,0) = Z653 E + 06(K~ 5 56 )9.62 E - 09(F 2 714 \0 + D^ 53,T 
= 2.552 E - 02(K -2 846 )(0 + 1)~° 5316 


Now, we can calculate the ISS Propulsion Module meteoroid PNP for an effective area of 
40.65 m 2 as: 


'90 79.8 


PNR 
= e 


40.65(^)J j2.552£-02(F- 2846 )(6>+l) 


- 0.5316 


dVdO 


= e 


0 18.8 


Meteoroid 

-4.4562T-04 = (^99955 


Note that the factor of 1/81 in front of the integral is necessary to guarantee that: 


. 90 79.8 

lj \f\VWde = 1 

O 1 n ley 


The equivalent meteoroid PNP using traditional risk assessment is 0.9993, resulting in a 
relative difference of +0.025%. 


49 





X ^ 

= s 
s 

2 E 
o . 

O S" 
0 ) - 2 : 

3 £ 

S o 

(0 

P 


8.00E-04 

7.00E-04 

6.00E-04 

5.00E-04 

4.00E-04 

3.00E-04 

2.00E-04 

1.00E-04 

0.00E-HD0 



-*• actual 
fractal 


0.000 0.200 0.400 0.600 0.800 1.000 


Critical Particle Diameter (cm) 


Figure 2.4.3-1. Local Meteoroid Flux Fractal for Relevant Whipple Shield 
Particle Diameters. Fractal Model is “Predicted.” ISS Propulsion Module 
Parameters (Meteoroid particles greater than 10 cm are tracked and 

avoided.) 
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Figure 2A.3-2. " t hstruc; Space" 2-D Fractal Representation of Whipple Shield 
Pertinent Meteoroid Flux for ISS Propulsion Module at Altitude of 394.7 km 
(Dimension = 3.333). Image created by replacing each line of a unit square by 10 lines, 
each 1/2 the length of the original, according to the formula r = (N) 1 D where r = scaling 
factor (2), N = number of lines replacing original line (10), and D = dimension of fractal 
(3.3333). Error = 0.25%. (Image created using BRAZIL DESIGN software.) 
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Figure 2.4.3-3. Fractal Model for Meteoroid Velocity Probability Distribution. 

ISS Propulsion Module Parameters 


2.4.4 Meteoroid Risk Assessment - Stuffed Whipple Shield Fractal Model 

In this section, a risk assessment is performed for the ISS Propulsion Module with Stuffed 
Whipple Shield configuration subject to the meteoroid. Figure 2.4.4- 1 shows the fit for 
the following “local” meteoroid flux fractal: 


F m (d) = 5.320E - 06(J~ 40097 ) 


This fractal was developed by considering the reduced range of values for meteoroid 
critical particle diameter pertinent to the Stuffed Whipple Shield configuration of the ISS 
Propulsion Module By considering only these particle diameters, a more accurate “local 
meteoroid flux fractal can be obtained for analysis purposes. The correlation (R) for this 
fractal is 1.0, with significance level of 0 0, indicating an excellent statistical fit. (Note that 
all R values reported are for the ln(Y) = ln(a) + b ln(X) model, not for the Y = a X 
model ) 

No graphical representation is provided for two reasons. The first is that this “fractal’ has 
almost integer dimension The second reason is that the dimension is much higher than a 
2-D representation could capture. 


Integration: 




We first combine the Stuffed Whipple Shield Phenomenology fractal 


D cru = 6.81F 


-0.5319 


(e + o' 


.26791 


with the local meteoroid flux fractal 


F m (d) = 5.320E- 06(J 40 ^ 7 ) 


to obtain the following (setting d - D cr jt): 


F s , (V,6) = 5.320 E - 06[6.81(< / ~ os: " 9 )(0 + 1)”"’" ] 
= 2.4281 £ - 09(K 2 1,5 )(<? + IT 1 0742 


0.2679 t -4.0097 


Next, we multiply this by the fractal model for meteoroid velocity probability distribution: 


f\V)F S{ (V,0) = 2.653 E + 06(K" 5 56 )2.4281E - 09 {V 2m )(6> + 1) 


,- 1 . 074 : 


= 6M2E-O3{V- iA21 )(0 + \) 


,- 1.0742 


Now, we can calculate the 1SS Propulsion Module meteoroid PNP for an effective area of 
40.65 m 2 as: 


90 79.8 


PNP = e 

1 1 ' 1 Meteoroid 

= e ' i(> = 0.9999960 


-40.65(— )f J6,442£-O3(r-’ J27 )(0+1) ' mn JVJ0 
^ 0 18.8 


Note that the factor of 1/8 1 in front of the integral is necessary to guarantee that: 


90 79.8 


81 


| \f(V)dVdO = \ 


0 18.8 


The equivalent meteoroid PNP using traditional risk assessment is 0.999982, resulting in a 
relative difference of +0.0014%. 
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Figure 2.4.4-1. Local Meteoroid Flux Fractal for Relevant Stuffed Whipple Shield 
Particle Diameters. Fractal Model is “Predicted." ISS Propulsion Module Parameieis 


2.4.5 Fractal Morphological Studies 

Figure 2.4. 5-1 shows a “morphing" transformation from the global orbital debris fractal 
(dimension = 2.1308) to the local orbital debris fractal (dimension = 1.4164) pertinent to 
the Whipple shield configuration of the ISS Propulsion Module. Given that this 
transformation occurs from a greater-than-2-D fractal to a lesser-than-2-D fractal, it seems 
reasonable that the global orbital debris fractal would need to lose some density. This 
density loss is evident in the second frame, and by the third and fourth frames, Ford [2001] 
observed that the mass making its way down from the center of the fractal is composed of 
images self-similar to the overall dynamic fractal. Furthermore, Ford has observed that the 
self-similarity continues “back into the larger portion of the fractal in a number of 
directions, and that the mass breaking off tends toward the lower right portion of the 
destination (local debris) fractal While this transformation is not entirely symmetrical in 
nature, these observations do lead one to conclude that the dynamics of morphing can be 
carried out with as much self-similarity as possible. 

Figure 2.4. 5-2 shows a morphing transformation from the global meteoroid fractal 
(dimension = 1.1782) to the phenomenology meteoroid fractal (dimension = 0.8143) for 
the Whipple Shield configuration of the ISS Propulsion Module. Because these 
dimensions are “closer" than those for the previous morphological study, it might be 
expected that less mass would need to be lost from the initial fractal to achieve the 
destination fractal. Again, it is evident by the third or fourth frame that this dynamic 
process is taking place in a self-similar fashion. 
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Figure 2.4.5-1. Transformation of Global Orbital Debris Fractal to Local (Whipple) 

Orbital Debris Fractal 
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Figure 2A.5-2. Transformation of Global Meteoroid Fractal to Phenomenology 

(Whipple) Meteoroid Fractal 
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2.4.6 “Finite Segment” Orbital Debris Risk Assessment - Stuffed Whipple Shield 

In this section, a “finite segment” representation of the ISS Propulsion Module has been 
approximated for the Stuffed Whipple Shield configuration subjected to orbital debris 
impacts. The analysis has been performed using an innovative spreadsheet approach to the 
traditional risk assessment, as well as an analogous fractal methodology 

The approach used is to segment the ISS Propulsion Module in an angular fashion around 
the circumference of the cylindrical portion Critical particle diameters are then calculated 
for each segment of the cylinder These critical particle diameters depend on one of two 
dominating impact angles. Figure 2.4.6- 1 provides a sketch of these- considerations. 

Figure 2 4.6-2 shows the “thermal” plot of the front and side views of the ISS Propulsion 
Module for the traditional risk assessment. The redder areas correspond to those areas or 
segments with lower relative PNP estimates, and thus, higher risk. The lighter yellow 
areas near the top and bottom of the Module reflect higher PNP estimates, and thus, lower 
risk. As expected, the front surface and angular segments near the center of the Module’s 
side have greater risk estimates. It is important to note that the radial symmetry of risk 
assessment is due to a lack of consideration of self-shadowing and shadowing by other 
ISS Elements. Thus, this “finite segment” approach is a first approximation to the ISS 
Propulsion Module risk assessment 

The overall PNP estimate using the “finite segment” approach is 0.9962, as compared to 
0.9965 for the “global” traditional risk assessment. 

Figure 2. 4. 6-3 shows the “thermal” plot of the front and side views of the ISS Propulsion 
Module for the fractal risk assessment To create this plot, fractal models were created 
for the RAM face, and each of the angular segments around the circumference of the 
cylindrical portion of the ISS Propulsion Module. These fractal models were then used to 
calculate PNP’s in a similar fashion as for the “global” assessments described in previous 
sections. 

In Figure 2.4 6-3, the redder areas again correspond to those areas or segments with 
lower relative PNP estimates, and thus, higher risk. The lighter yellow areas near the top 
and bottom of the Module reflect higher PNP estimates, and thus, lower risk As 
expected, the front surface and angular segments near the center of the Module’s side 
have greater risk estimates It is important to note that the radial symmetry of risk 
assessment is due to a lack of consideration of self-shadowing and shadowing by other 
ISS Elements Thus, this “finite segment” approach is a first approximation to the ISS 
Propulsion Module risk assessment. 

The overall PNP estimate using the “finite segment” approach is 0.9955, identical to that 
for the “global” fractal risk assessment. 
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Figure 2.4.6- 1. Sketch of “Finite Segment” Approach 
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Figure 2.4.6-2. “Finite Segment” Results - Traditional Assessment (Stuffed 
Whipple Shield Configuration - Orbital Debris) 
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Figure 2.4.6-3. “Finite Segment” Results - Fractal Assessment (Stuffed Whipple 

Shield Configuration - Orbital Debris) 


3.0 Conclusions 


1 The meteoroid and orbital debris environment definitions pertinent to the ISS 
Propulsion Module fit fractal forms with a high degree of visual and statistical 
confidence 

2. The ORDEM96 orbital debris environment underestimates the flux relative to the 
NASA TM-4527 orbital debris environment definition for the ISS Propulsion Module 

3. Rectangular cross-sectional projected areas are self-similar and can be reduced and 
dealt with accordingly While not strictly a fractal form, the self-similarity of these 
areas can be dealt with in a manner consistent with fractal forms 

4 Ballistic limit and crater phenomenology models pertinent to the ISS Propulsion 
Module fit fractal forms with a high degree of visual and statistical confidence. 
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5. Fractal risk assessments of ISS Propulsion Module orbital debris PNP for Whipple and 
Stuffed Whipple shielding configurations very closely estimate PNP values obtained by 
traditional spreadsheet analysis. 

6. Fractal for Whipple Shield phenomenology model is excellent fit with correlation of 
0.991 and significance level of 4 x 10 54 

7. Fractal for “local” orbital debris flux model is excellent statistical and visual fit with 
correlation of 0.967 and significance level of 5.6 x 10 6 

8. Non-fractal approximating model for Whipple Shield orbital debris velocity probability 
distribution is excellent statistical fit with correlation of 0.985 and significance level of 
1.3 x 10* 45 . Despite excellent statistical fit, non-fractal approximating model does not 
display good visual fit. 

9. Revised preliminary standard (non-ffactal) risk assessment indicates the probability of 
no orbital debris penetration (PNP) for ISS Propulsion Module Whipple Shield 
configuration is approximately 0 9876. 

10 Preliminary fractal risk assessment for ISS Propulsion Module Whipple Shield 
configuration results in orbital debris PNP of 0.9891, a difference of 0. 1 5 /o from the 
equivalent standard assessment. 

1 1 Fractal for Stuffed Whipple Shield phenomenology model provides excellent statistical 
and visual fits with correlation of 0.9766 and significance level of 9.5 x 10' 4 “. 

12. Fractal for “local” orbital debris flux model for Stuffed Whipple Shield is excellent 
statistical and visual fit with correlation of 0.9307 and significance level of 6 8 x 10 

13. Revised preliminary standard (non-fractal) risk assessment indicates the probability of 
no orbital debris penetration (PNP) for ISS Propulsion Module Stuffed Whipple 
Shield configuration is approximately 0.9965. 

14. Preliminary fractal risk assessment for ISS Propulsion Module Stuffed Whipple Shield 
configuration results in orbital debris PNP of 0.9955, a difference of 0.10% from the 
equivalent standard assessment. 

15. Fractal crater phenomenology model shows excellent statistical fit with correlation of 
0.922 and significance level of 3.03E-09. 

16. Fractal and traditional spreadsheet methodologies can result in meaningful “finite 
segment” risk assessments. Differences between the “finite segment” and “global risk 
assessment PNP values are insignificant in both cases. 
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4.0 Recommendations 


1 Use the TM-4527 environment model definitions with constant orbital debris particle 
density of 2.8 gm/cnr for Propulsion Module requirements and consistency with 
BUMPER 

2. Protect the ISS Propulsion Module to the same extent as the manned modules by 
using the Stuffed Whipple Shield design. 

3. Perform fractal sensitivity analyses for other geometries, orbit parameters, and shield 
configurations 

4. Develop a generic fractal M/OD analysis spreadsheet 
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